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Abstract 

The running quark masses mq(/i) at various energy scales /i (p = 
IGeV, p = rUg, fi = mz and so on) are evaluated by using the mass 
renormalization equations systematically. Also, those at energy scales fi 
higher than fi = mz (from fi = 10^ GeV to fi = 10^® GeV) are evaluated 
by using the evolution equations of the Yukawa coupling constants for 
the standard model with a single Higgs boson. 
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1. Introduction 


It is very important to know reliable values of quark masses rUq not only for 
hadron physicists but also for quark-lepton physicists. For such a purpose, a review 
article by Gasser and Leutwyler [1] has offered useful information on the running 
quark masses mq(/i). However, since Gasser and Leutwyler’s review article [1] in 
1982, there have been some developments, for example, higher order calculation of 
perturbative QGD [2, 3], matching condition at quark threshold [4], and discovery 
of the top quark [5]. 

In this talk we report the running quark masses at various energy 

scales jj, {jj, = IGeV, ^ = niq, jj, = mz and so on) which are evaluated by using 
the mass renormalization equations systematically. The calculation was done by 
taking the matching condition at the quark flavor threshold into account. Also, 
those at energy scales fx higher than n = mz (from jj, = 10^ GeV to /i = 10^® GeV) 
are evaluated by using the evolution equations of the Yukawa coupling constants 
for the standard model with a single Higgs boson. 

In the next section, we review values of light quark masses mu(/i), 
and ms(/i) at /i = 1 GeV. In Sec.3, we review pole mass values of heavy quark 
masses and In order to estimate mq{^) at any /i, we must 

know the values of the QGD parameters (n = 3,4, 5, 6). In Sec.4, the values of 
Q;s(/i) and A^ are reviewed. In Sec.5, running quark masses mq{fi) are evaluated 
for various energy scales fx, e.g., fx = 1 GeV, fx = mq, fx = fx = mz, 

fx = Aw, and so on, where is a “pole” mass of the quark, and A^ is a 

symmetry breaking energy scale of the electroweak gauge symmetry SU(2)£, x 
U(l)y, i.e. Aw = (0°) = {^\/2Gf)~^/^/2 = 174.1GeV. In Sec.6, the reliability of 
the perturbative calculation below /i ~ 1 GeV is discussed. In Sec.7, evolution of 
the Yukawa coupling constants of the standard model with a single Higgs boson 
is estimated for energy scales higher than g, = Aw- Finally, Sec.8 is devoted to 
summary and discussion. 

2. Review: light quark masses at /i = 1 GeV 

Since Gasser and Leutwyler [1] have obtained the light quark masses m„(/i), 
md{fx) and ms{fx) at fx = 1 GeV, various values of light quark masses are reported. 
We summarize these values in Table 1. 

As shown in Table 1, there is not so a large discrepancy among these esti¬ 
mates as far as and md are concerned. But, for the strange quark mass ms, 
two different values, mg — 175 MeV and mg — 200 MeV have been reported. We 
use weighted averages as input values in our calculations. 
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Table 1 . Light quark mass values at 1 GeV (in unit of MeV) 



mu 

ma 

ms 

Gasser and Lentwyler (1982) [1] 

5.1 

±1.5 

8.9 

±2.6 

175 

±55 

Dominguez and Rafael (1987) [6] 

5.6 

±1.1 

9.9 

±1.1 

199 

±33 

Narison (1995) [7] 

4 

±1 

10 

±1 

197 

±29 

Leutwyler (1996) [8] 

5.1 

±0.9 

9.3 

±1.4 

175 

±25 

Weighted averages 

4.90 

±0.53 

9.76 

±0.63 

187 

±16 


3. Review: pole masses of heavy quarks 

Charm and bottom quark masses 

Gasser and Leutwyler (1982) [1] have estimated charm and bottom quark 
masses rric and rrih and Tirard and Yndnrain (1994) [9] have re-estimated rric and 
mb precisely and rigorously. On the other hand, from ip- and T-sum rules, Narison 
(1994) [10] has estimated the running qnark masses corresponding to the short- 
distance pertnrbative pole masses to two-loops and three loops. In Table 2, we 
summarize their values. We use weighted averages in Table 2 as inpnt valnes in 
onr calculations. 


Table 2 . Pole masses of charm and bottom quark 



MPole 

MPole 

Tirard and Yndnrain (1994) [9] 

1.570 

±0.019^0.007 

4.906 

+0.069 -p ri riri/i+O.Oll 

-0.051 + U.UUU_Q Q4Q 

Narison (1994) [10] 

1.64 

± 0.03 

4.87 

±0.05 ±0.02 

Weighted averages 

1.59 

±0.02 

4.89 

±0.05 


Top quark mass 


The top quark mass values obtained by 
by the CDF collaboration (1994) [5, 11] 
and the DO collaboration [12] are sum¬ 
marized in Table 3. We use the val¬ 
ues quoted by the particle data group 
(PDG96) [13] as the pole mass of the top 
quark. 


Table 3 . Pole mass of top quark 



MPole 

CDF (1994) [5] 

174 ±10 

CDF (1995) [11] 

176 ±8 ±10 

DO (1995) [12] 

199 ±22 

PDG (1996) [13] 

180 ±12 
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4. Estimation of as{i^) and 

Prior to estimates of the running quark masses we must estimate the 

values of Q;s(/i) and A^. The effective QCD coupling o;* = fi'^/dTr is given by [14] 


(h) 


47r 1 

[ 2(3i\nL 

, 4/32 


ToL' 

1 (^l L 

AV2 

_V 2 ) 8 / 3 ? 4_ 





where 


n/ \ /^O 2 A 

/3(».) = 


3 <^2 4 


a, — 


6471-3 


+ 0(q;^) , 


2 

-1 

3 


19 

y 


5033 


325 


/^o = 11 - 77^g , A = 51 - —Uq , (32 = 2857- —Uq + — 


(4.1) 

(4.2) 

(4.3) 


L = ln(/iVA2) . (4.4) 

The value of ag^jj) is not continuous at nth quark threshold /i„ at which the 
nth quark ffavor channel is opened, because the coefficients /3o, A and P 2 depend 
on the effective quark ffavor number Uq. 


Therefore, we use the ex¬ 
pression Q;("')(/i) with a dif¬ 
ferent A^ for each energy 
scale range /in < /W < f^n+i- 
The values of A^ are evalu¬ 
ated by matching condition 
[4]. In Table 4, the val¬ 
ues of A^ are summarized 
and the underlined values 
denote input values A^. 
We show the threshold be¬ 
havior of Q;^"'^(/i) in Fig. 1. 
We can see that 
in /i„_i < /i < /in connects 
with in /in < /i < 

/in+i continuously. 



Fig. 1. The threshold behavior of 
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Table 4 . The values of in unit of GeV. 


n 

3 

4 

5 

6 

A'ls 

n qqq+0.047 
U.OOO_0.042 

0.2911^^? 

0.209i:]"^^ 

0.0882tHli 


5. Estimation of running quark masses mq{mq) 

From the pole mass values we estimate the running mass values mq{n) 

at /i = for q = u,d, s hj using the relation [15] 


rriq (Mf= Mf 


1 + 


3 TT V TT / 


+ 


(5.1) 


The values of K, and mq^jj) at for q = c,b,t are summarized in Table 5. 

Table 5 . Parameter K , Pole masses and Running mass mq{fi) at 



K 

MPole (QgY) 

mq{MP°^^) (GeV) 

c 

14.47 

1.59 

1.213 

y 

12.94 

4.89 

4.248 

t 

10.98 

180 

170.1 


The scale dependence of a running quark mass mq{n) is governed by the 
equation [2] 


A 

djj, 


IJ—m,(p) = -l(a,)7n,(p) , 


(5.2) 


where 


7(«.)=7o- + 7i(-)' + 72(-)' + 0(«^) , (5.3) 

71 \ 71 / \ 71 / 


101 5 1 

70 = 2 . = .72 = ^^ 


/2216 160,, A 

- {-w + —‘^<=*4"« 


140 2 ' 


The running quark mass mq{fi) is given by 

rriqiiJ,) = R{as)mq , 


(5.4) 

(5.5) 


R{as) = 


/3i7o\ 


2 TT 


1 + 2 


a, 


Po Pi 


5 





2 


(5.6) 


+ 


1 

2 



Pao \ 

Pi ) 



Pi7i 

Pi 


P2l0 PllO 
16/32 ^ 2/33 



where niq is the renor¬ 
malization group invariant 
mass which is independent 
of ln(/r2/A2), tts is given by 
(4.1) and {3i {i = 0,1, 2) are 
defined by (4.3). 

By using A^, we can eval- 

iVlo 

uate i?*'"^(/i) for /i < 
where fin is the nth quark 
flavor threshold and we take 
= mqnimqn). We shoW 
the threshold behavior of 
R^'^Pfi) in Fig. 2. As shown 
in Fig. 2, the behavior of 
R{fi) is discontinuous at 
h /^n = '^qni^q-n)- 



Fig. 2. The threshold behaviour of 


We can evaluate the values of mq{mq) (g = c,b,t) by using the values of 
and the relation 


^5n(/^) 




(5.7) 


Similarly, we evaluate the light quark masses uiqiniq) by using the values 
mg(lGeV) and the relation 


= i?*'^^(/i)/i?^^^(lGeV) mg(lGeV) (/i < p.4) 


(5.8) 


Running quark mass values rriqnifi) at /i > fin+i cannot be evaluated by using 
R{fiY'^^ straightforwardly, because of the discontinuity of R{fi) at quark threshold 

fi = fln = mq^{mqP). 

The behavior of the nth quark mass {n < N) aX fi^ < f^ < f^N+i is 
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given by the matching condition [16] 



‘I#*)/ 




5 89\ 



(5.9) 


where 

a;AT = In . (5.10) 

In Fig. 3, we illustrate the /i-dependency of the light quark masses rriq^jj) (g = 
M, d, s) which take the matching condition (5.9) into account. We also illustrate 
the behavior of the heavy quark masses mg(/i) (g = c, b, t) in Fig. 4. 



H(GeV) 


^(GeV) 


Fig. 3. running masses of Fig. 4. running masses of 

light quarks heavy quarks 

The numerical results are summarized in Table 5, where input values rriq (1 
GeV) for q = u,d,s and for q = c,b,t are used. The hrst and second 

errors come from (or and ±AA^, respectively. The values with 

asterisk should not be taken rigidly, because these values have been calculated in 
the region with a large 

6. Reliability of the perturbative calculation below^ /i ~ 1 GeV 

As we noted already, the values of the light quark masses mq{mq) (g = m, d, s) 
should not be taken rigidly, because the perturbative calculation below /i ~ 1 GeV 
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Table 5 . Running quark mass values rriq^jj) at /i = fUq 


q = 

u 

d 

S 

C 

b 

t 

MPole 

*0.501 

*0.517 

*0.681 

1.59 

4.89 

180 


±0.002 

±0.002 

±0.011 

±0.02 

±0.05 

±12 


+ 0.058 

- 0.052 

+ 0.067 

- 0.060 

+ 0.065 

- 0.059 




mq{MP°^^) 

*0.0308 

*0.0443 

*0.275 

1.213 

4.248 

170.1 

+ 0.0018 

- 0.0019 

+ 0.0016 

- 0.0017 

+ 0.015 

- 0.016 

±0.018 

±0.046 

±11.4 


+ 0.0002 

- 0.0005 

+ 0.0001 

- 0.0006 

- 0.004 

+ 0.001 

- 0.040 

+ 0.034 

- 0.040 

+ 0.036 

^0.3 

mq{mq) 

*0.436 

*0.448 

*0.549 

1.302 

4.339 

170.8 


+ 0.001 

- 0.002 

±0.001 

±0.007 

±0.018 

±0.046 

±11.5 


+ 0.058 

- 0.052 

+ 0.059 

- 0.053 

+ 0.059 

- 0.052 

- 0.020 

+ 0.019 

- 0.029 

+ 0.027 

^0.2 

mq{lGeV) 

0.00490 

0.00976 

0.187 

1.467 

6.356 

339 


±0.00053 

±0.00063 

±0.016 

±0.028 

±0.080 

±24 





+ 0.005 

- 0.002 

+ 0.214 

- 0.164 

+12 

-11 

mq{mc) 

0.00421 

0.00838 

0.160 

1.302 

5.782 

318 

rric = 1.302 

±0.00045 

±0.00054 

±0.014 

±0.018 

±0.047 

±22 


- 0.00011 

+ 0.00007 

- 0.00021 

+ 0.00015 

- 0.004 

+ 0.003 

- 0.020 

+ 0.019 

+ 0.145 

- 0.112 

+10 

-9 

mq{mb) 

0.00312 

0.00621 

0.119 

0.950 

4.339 

253 

rrib = 4.339 

±0.00034 

±0.00040 

±0.010 

±0.016 

±0.046 

±18 


- 0.00020 

+ 0.00016 

- 0.00040 

+ 0.00031 

- 0.008 

+ 0.006 

- 0.052 

+ 0.045 

- 0.029 

+ 0.027 

+4 

-3 

mq{mw) 

0.00224 

0.00446 

0.0855 

0.668 

3.029 

182 

ruw = 80.33 

±0.00024 

±0.00029 

±0.0073 

±0.013 

±0.038 

±13 


- 0.00017 

+ 0.00014 

- 0.00035 

+ 0.00029 

- 0.0066 

+ 0.0055 

- 0.047 

+ 0.043 

- 0.074 

+ 0.069 

±0.04 

mq{mz) 

0.00222 

0.00442 

0.0847 

0.661 

2.996 

180 

mz = 91.187 

±0.00024 

±0.00029 

±0.0072 

±0.012 

±0.038 

±13 


- 0.00017 

+ 0.00014 

- 0.00034 

+ 0.00029 

- 0.0066 

+ 0.0055 

- 0.047 

+ 0.042 

- 0.074 

+ 0.069 

±0.02 

mq{nit) 

0.00212 

0.00422 

0.0809 

0.630 

2.847 

170.8 

rrit = 170.8 

±0.00023 

±0.00027 

±0.0069 

±0.009 

+ 0.021 

- 0.020 

±11.5 


- 0.00017 

+ 0.00014 

- 0.00033 

+ 0.00028 

- 0.0063 

+ 0.0053 

- 0.045 

+ 0.041 

- 0.074 

+ 0.069 

^0.2 

mq{Aw) 

0.00212 

0.00422 

0.0808 

0.629 

2.843 

170.5 

Aw = 174.1 

±0.00023 

±0.00027 

±0.0069 

±0.012 

±0.036 

±12.3 


- 0.00017 

+ 0.00014 

- 0.00033 

+ 0.00028 

- 0.0063 

+ 0.0053 

- 0.045 

+ 0.041 

- 0.075 

+ 0.070 

^0.3 


















































































seems to be not very reliable. Let us look at this more explicitly. 

In order to see the reliability of the calculation of as{fi) by using (4.1), in 
Fig. 5, we illustrate the values of the second and third terms in { } of (4.1) 

separately. The values of the second and third terms exceed one at /x ~ 0.5 GeV 
and n ~ 0.6 GeV, respectively. Also, in Fig. 6, we illustrate the values of the 
second and third terms in { } of (5.6) separately. The values of the second and 
third terms exceed one at p ~ 0.6 GeV and p ~ 0.7 GeV, respectively. These mean 
that the perturbative calculation is not reliable below p ~ 0.7 GeV. Therefore, the 
values with asterisk in Tables 5 should not be taken strictly. 




Fig. 5. the values of the sec- Fig. 6. The values of the sec¬ 
ond and third term of (4.1) ond and third term of (5.6) 

7. Evolution of Yukawa coupling constants 

By using the renormalization group equation, we estimate the Yukawa cou¬ 
pling constants in the standard model with a single Higgs boson. 

The quark mass matrices and at /i = Ky/ are given by 

^a(h) = , (7.1) 

where Ya denotes a matrix of the Yukawa coupling constants yfj (a = u,d-,i = 
1, 2, 3), (Ya)ij = ylj and v is the vacuum expectation value of the Higgs boson. 
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The renormalization scale dependence of a matrix Ha = YaY^ is given by [3] 


dt 


Ha = 






(IGtt 






Ha + H, 


IGtt^ 


;3<‘>' + 


(IGtt 


2'|2^a 




(7.2) 


where t is given by t = \n{n/mz) and the one-loop and two-loop contributions 
(3^^ and g^Ye written as (3^'^ = -|- and (3j3) = Yjb^\Hb + 

Y.b,cb'^HbHc, respectively. The expressions of coefficients cb\ a^, etc. have been 
given in Ref. [3]. For the input parameters, we use the quark masses mq{mz) in 
Table 5 and the following parameters in the CKM matrix V [13]: 

\Vus\= 0.2205 ± 0.0018, \Vcb\ = 0.041 ± 0.003, iKfe/Kbl = 0.08 ± 0.02 . (7.3) 
For the gauge coupling constants, we use [17] 


a{mz) = (128.89 ± 0.09)-^ sin^ 9w = 0.23165 ± 0.000024, 
asimz) = 0.118 ±0.003 . 


(7.4) 


The input value of Higgs boson uih is V^Aw = 246.2 GeV. 


If the input value mn is less 
than 2.2 x lO^ GeV (2.3 x 
10^ GeV) for two (one) loop 
evaluations, then the quar- 
tic coupling constant A, of 
the Higgs boson self inter¬ 
action, becomes negative at 
high energy. On the other 
hand, if the input value 
rriH is more than 2.6 x 10^ 
GeV for both two and one 
loop calculations, the burst 
of A occurs at high en¬ 
ergy. In Fig. 7, we illustrate 
the p-dependency of the 
Yukawa coupling constants 
Vqifi) {q = u,d,s,c,b,t) 
which take the renormaliza¬ 
tion equation into account. 



u 


10 ° 10 "^ 10 ° 10 ^^ 10 ^° 
p(GeV) 


Fig. 7. Evolution of the Yukawa cou¬ 
pling constants 
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8. Summary 

We have evaluated the running quark mass values TUq^jj) at various energy 
scales below /x = Aw and the Yukawa coupling constants of the standard model 
with a single Higgs boson at energy scales above fx = Aw- Although we have used 
the renormalization equation, the perturbative calculation below /i ~ 0.7 GeV is 
not adequately reliable because the values of the second and third terms in the { 
} of perturbative series (4.1) and (5.6) exceed one less than fx ~ 0.7 GeV. 

We discuss the grade of parameters htted in mass matrix models. At present, 
the “conhdence” grade of the “observed” values of the running quark masses mq^jx) 
and GKM matrix parameters are not equal at the same levels because these val¬ 
ues are highly dependent on models or other experimental values (input values). 
Therefore, it is important for the model-building of quark mass matrix that we 
know the conhdence levels of these values. Our opinion based on the present work 
is summarized in the following Table: 


Table 6 . The reliability of the “observed” values 


grade 

GKM matrix element 

quark mass ratio 

(i) 

(Reliable) 

IW.I 


(ii) 

(Almost reliable) 

iWftI 

md/ms mdmb mb/mt 

(hi) 

(Somewhat variable) 

iKfel 

mc/mt mu/me ms/me 

(iv) 

(Variable) 


mu/md md/mt mg/mb 

(v) 

(Unreliable) 

\Vtd\ 



We have classihed the GKM matrix elements on the basis of the experimental and 
theoretical errors. In grading the quark mass ratios, we have considered ratios are 
reliable in the cases where (1) both two quarks are heavy quarks or light quarks 
and (2) the mass difference between two quarks is small. 

Finally, we would like to point out that we should use the running mass 
values of /i = mz rather than /i = 1 GeV for quark mass matrix phenomenology, 
together with the GKM matrix parameters at /i = mz- 
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